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The origin of the 90 K nematic transition in the chalcogenide FeSe, which displays no magnetic 
order down to T = 0, remains a major puzzle for a unifying theory for the iron-based superconduc¬ 
tors. We analyze this problem in light of recent experimental data which reveal very small Fermi 
pockets in this material. We show that the smallness of the Fermi energy leads to a near-degeneracy 
between magnetic fluctuations and fluctuations in the charge-current density-wave channel. While 
the two fluctuation modes cooperate to promote the same preemptive Ising-nematic order, they 
compete for primary order. We argue that this explains why in FeSe the nematic order emerges 
when the magnetic correlation length is smaller than in other Fe-based materials, and why no mag¬ 
netism is observed. We discuss how pressure lifts this near-degeneracy, resulting in a non-monotonic 
dependence of the nematic transition with pressure, in agreement with experiments. 


Nematic order in Fe-pnictides and Fe-chalcogenides de¬ 
velops at a temperature Tg that is larger than the mag¬ 
netic transition (for reviews, see Q)- It spontaneously 
breaks the tetragonal C 4 lattice symmetry down to or¬ 
thorhombic C2. The origin of this symmetry breaking 
is currently one of the most intensely debated issues of 
the Fe-based superconducting materials 0. In the Fe- 
pnictides, nematic order occurs reasonably close to the 
instability towards stripe magnetic order at the Neel tem¬ 
perature Tjv. Because the stripe order breaks Z 2 tetrag¬ 
onal symmetry {C 4 —iC 2 ) in addition to the 0(3) spin- 
rotational symmetry and because Tg and Tjv show similar 
doping dependencies, it seems reasonable to associate the 
nematic order with magnetism Q . Indeed, several groups 
have argued that magnetic fluctuations split the 

mean-field stripe magnetic transition into two separate 
0(3) and Z 2 transitions. The discrete Z 2 symmetry is 
broken first at Tg > Tjv, resulting in an intermediate 
phase, dubbed Ising-nematic, where long-range magnetic 
order is absent but the O4 lattice symmetry is broken 
down to O2. Such Z 2 order triggers orbital and struc¬ 
tural order as all three break the same O4 symmetry. 


The magnetic scenario for nematicity in Fe-pnictides 
is supported by a variety of experimental observations, 
such as the doping dependencies of T^r and Tg [l^, the 
scaling between the shear modulus and the spin-lattice 
relaxation rate [l^, and the sign-change of the in-plane 
resistivity anisotropy between electron-doped and hole- 
doped Fe-pnictides [ij. This scenario, however, has been 
challenged for the Fe-chalcogenide FeSe. This material 
displays a nematic transition at Tg Ki 90K. The prop¬ 
erties of the nematic phase in FeSe resemble those in 
Fe-pnictides: similar softening of the shear modulus M, 
similar orthorhombic distortion and orbital order iB-flsj. 
and similar behavior of the resistivity anisotropy upon 
applied strain [T^ . Furthermore, neutron scattering ex¬ 
periment shows that spin fluctuations are peaked at the 
same ordering vectors as in the Fe-pnictides [23,[ll|. Yet, 
in distinction to Fe-pnictides, no magnetic order has thus 
far been observed in FeSe in the absence of external pres¬ 


sure [2^ [2^. Moreover, NMR measurements were inter¬ 
preted as evidence that the magnetic correlation length 
^ remains small at Tg [l^, [2^. Although in the Ising- 
nematic scenario ^ does not have to be large at Tg, this 
seems to be the case for all Fe-pnictides. 

Given these difficulties with the Ising-nematic sce¬ 
nario, spontaneous orbital order has been invoked to ex¬ 
plain the nematic state in FeSe mill. However, at 
present, no microscopic theory exists where orbital order 
appears spontaneously instead of being induced by mag¬ 
netism |25l - l^ . Alternative scenarios for magnetically- 
driven nematicity in FeSe have also been proposed, in¬ 
volving the formation of a quantum paramagnet (30j| , the 
onset of spin quadrupolar order 13111 . and strong frustra¬ 
tion of the magnetic fluctuations [^. Yet, the issue of 
why FeSe does not fit into a “universal” theory for the 
iron-based superconductors still persists. 

In this communication, we present an extension of the 
spin-nematic scenario which explicitly builds on a unique 
property of the electronic structure of FeSe, namely, the 
fact that the Fermi energy Ep in this material is small 
- only a few meV, as seen by ARPES and dH-vA ex¬ 
periments [Hill. For a system with a small Ep, ear¬ 
lier renormalization-group (RG) calculations have shown 
that there are two density-wave channels whose fluctua¬ 
tions are strong at momenta (0, 0): a spin density- 

wave (SDW) channel and a charge-current density-wave 
(CDW) channel (a CDW with ima gina ry order param¬ 
eter, which we denote as iGDW [3^). The relative 
strength between the two depends on the sign of the inter¬ 
pocket exchange interaction (u 2 in our notations below). 
For repulsive U 2 , the coupling in the SDW channel is 
larger, while for attractive U 2 the coupling in the iCDW 
channel is larger. In both cases, however, the RG calcu¬ 
lations show that the coupling in the subleading channel 
approaches the one in the leading channel at small ener¬ 
gies. The RG process stops at Ep, implying that if Ep 
is larger than the highest instability temperature (Tg for 
FeSe) the subleading channel is not a strong competitor 
and for all practical purposes can be neglected. However, 
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FIG. 1: Schematic representation of the SDW and iCDW 
ordered states with ordering vector (7r,0). In the latter, the 
arrows represent charge currents along the bonds, while in the 
former, they represent spins on the sites. While fluctuations 
of both channels support nematicity, they compete for long- 
range magnetic and charge order. 


if Ep ~ Tg, as in FeSe, the couplings in the two chan¬ 
nels become degenerate within the RG. The degeneracy 
implies that the order parameter manifold increases from 
0(3) X Z2, for the three-component SDW, or from Z2 x Z2, 
for the one-component iCDW, to a larger 0(4) x Z 2 . In 
all cases, the Z 2 part of the manifold corresponds to se¬ 
lecting either (0 , tt) or (tt, 0 ) for the density-wave ordering 
vector. While in both 0(3) x Z 2 and 0(4) x Z 2 models 
the Z 2 symmetry can be broken before the continuous 
one, in the latter this happens at a significantly smaller 
correlation length. As a result, at small Ep, the nematic 
order emerges while magnetic fluctuations are still weak. 
Furthermore, the SDW transition temperature in the 
0(4) model is additionally suppressed due to the compe¬ 
tition with iCDW. We argue that these features explain 
the properties of the nematic state in FeSe, including 
non-monotonic pressure dependence of Tg [35|, l36| . 

The model. We consider a quasi-2D itinerant band 
model with two hole pockets at the F point and two elec¬ 
tron pockets at (0, tt ) and (tt, 0) in the 1-Fe Brillouin zone 
[lOL . This model can be obtained from an underly¬ 
ing 5-orbital model with Hubbard and Hund interactions 
and hopping between the Fe 3d orbitals 

The quadratic part of the Hamiltonian in the band ba¬ 
sis describes the dispersion of the low-energy fermions, 
and the information about the orbital content along the 
Fermi pockets is passed onto inter-pocket and intra¬ 
pocket interactions, which are the Hubbard and Hund 
terms dressed by the matrix elements ari sing from the 
change from the orbital to the band basis [4^. The an¬ 
gular dependence of the matrix elements leads to angle- 
dependent interactions. The three interactions rele¬ 
vant for Ising-nematic order are the inter-pocket density- 
density interaction ui, the exchange interaction U 2 , and 
the pair-hopping interaction U 3 14III. To simplify the 
analysis, we follow earlier works [l0| and analyze the 
Ising-nematic order within an RG procedure that (i) ap¬ 
proximates these three interactions as angle-independent 
and (ii) restricts the analysis to one hole pocket. The 
extension to two pockets and angle-dependent interac¬ 
tions makes the calculations more involved but does not 


modify the RG equations in any substantial way and, 
moreover, leaves them intact if we treat the hole pockets 
as circular and neglect the dxy orbital component on the 
electron pockets. 

We label the fermions near the hole pocket as Ck and 
the fermions near the electron pockets as /i,k and / 2 ,k. 
The 0(3) magnetic order parameter is given by 

= (4,a'^«/3/j,k+Q,,/3 + /j-c) , (I) 

whereas the Z 2 iGDW order parameter is 

~ ]V ~ h-c^ , (2) 

ka 

with j = 1,2 corresponding to the two possible ordering 
vectors Qi = ( 7 r, 0 ) and Q 2 = ( 0 , 7 r). We show these two 
ordered states in Fig. [T] 

0(4) Ising-nematic action. In the Ising-nematic sce¬ 
nario, the O 4 —>■ O 2 symmetry breaking implies the ap¬ 
pearance of a composite order, quadratic in the density- 
wave order parameters Mj and 4)j. To analyze this sce¬ 
nario, we need to know the flow of the couplings that 
drive SDW order, Fgdw = ui U 3 , and iCDW order, 
Ticdw — “1“ ^3 2 zi 2 (Ref. 0 ). The bare coupling 

Tsdw ^ bjcdw when 112 ^ 0 and Fjcdw ^ Tgdw when 
U 2 < 0. As one integrates out the high-energy degrees of 
freedom via an RG procedure, the ratio M 2 /(mi + M 3 ) de¬ 
creases as the system flows to lower energies (or temper¬ 
atures) and approaches zero at the energy/temperature 
scale in which the system develops SDW/iCDW order. 
This holds, however, only if this scale is larger than Ep. 
If Ep is larger, the RG flow stops at Ep and the sys¬ 
tem develops an instability only in the channel with the 
largest bare coupling. 

To illustrate our point, we plot in Fig. [3(a)-(b) the 
RG flow of Ficdw and Fgdw for a particular set of bare 
couplings Ml ( 0 ) = M 2 ( 0 ) = IOM 3 ( 0 ), chosen deliberately 
to give a negative bare Ficdw Under the RG flow, Fjcdw 
becomes positive and approaches Fgdw at the scale where 
the couplings diverge and the system develops a density- 
wave order. The Fermi energy Ep sets the scale at which 
the RG flow stops. In case I (large Ep), the RG stops 
when Ficdw is still small. In case H (smaller Ep), the RG 
stop when r^dw is comparable to Fgdw, and in case HI 
(even smaller Ep), the RG flow reaches the 0(4) fixed 
point already at energies larger than Ep. We associate 
case I in Fig. [5] with Fe-pnictides, and cases H/HI with 
FeSe based on the values ot Ep obtained by ARPES and 
quantum oscillations [mu. 

We next take the RG results as input and analyze 
the emergence of a nematic order which spontaneously 
breaks the symmetry between momenta Qi and Q 2 with¬ 
out breaking any other symmetry. The analysis follows 
the same steps as for pure SDW order [l^ : we introduce 
Mj and 4)^ {j = 1,2) as Hubbard-Stratonovich fields 
which decouple the four-fermion interaction terms, inte¬ 
grate over the fermions, and obtain the effective action 
in terms of Mj and : 
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FIG. 2: (a) RG flow of the SDW and iCDW interactions Fgdw 
(red curve) and Ficdw (blue curve) as function of decreasing 
energy E. W is the bandwidth, uq = ui (0) = U 2 (0) = 
10 u 3 (0) is the bare interaction parameter, and the dashed line 
is the energy in which the two degenerate instabilities occur. 
The RG flow stops at the Fermi energy Ep: if Ep is large (case 
I, Fe-pnictides), only SDW fluctuations are relevant, whereas 
if Ep is small (cases II/III, FeSe), both SDW and iCDW 
fluctuations are important. The insets show schematically 
the Fermi pockets in each case, (b) Ratio Ficdw/Fsdw along 
the RG flow, (c) Electronic manifestation of the Ising-nematic 
order on the hole pockets. There is a cos 26 distortion, with 
opposite signs for the two pockets, and an overall shift of the 
chemical potential. 


J qj ^ J xj 

- + + (3) 

where x7,l = KcL “ and x7,l = - n, with 

n, = /,G, ^k+q {Gf^^k + Gf 2 ,k)- Note, the only asymme¬ 
try between the two order parameters is due to the inter¬ 
actions Fsdw and Ficdw, respectively. Near Qj, we can ex¬ 
pand X7(l),q ~ Fo,^(c)-ha(q-Qj)^, where ro,^(c) measures 
the distance to the SDW (iCDW) mean-field instability 
and a ^ G>{1). The input from the RG analysis is that 
ro,s and rg^c are close to each other. The quartic coeffi¬ 
cients are given by (rt, g) = ±5 (G/^^fe ± Gf^^kf. 

At Fsdw = Ficdw, the action depends on M and <i> only 
via the combination -|- and the order parame¬ 
ter manifold is 0(4) x Z 2 . Evaluating the integrals at 
Ep ~ Tg, we find u > 0 and g > Q, what implies that 
long-range order selects either j = 1 or j = 2 , but not 
both, i.e. it breaks both 0(4) and Z 2 symmetries. 


Within a mean-field approximation, 0(4) and Z 2 are 
broken at the same temperature. Beyond mean-field, 
the Z 2 symmetry is broken first, and both M and <i> 
contribute to it, even if Fsdw ^ Ficdw- To see this, 
we treat M and $ as fluctuating fields, introduce the 
composite fields ■0 = u (M^ + -I- -f and ip = 

g (M^ + — M.y + to decouple the quartic terms, 

integrate over the primary fields M and $ and obtain the 
action in terms of ip and p: 


Ses [p, Ip] = 


p^ 02 3 r 

1 
2 


/ 


In 


(x0^+0) 


(4) 


The field ip has a non-zero expectation value {ip) yf 0 
at any tenperature as it does not break any symmetry, 
but only renormalizes the correlation lengths of the pri¬ 
mary fields M and 4) to = fo,s(c) + (0)- A non-zero 
((/?), on the other hand, breaks the tetragonal G 4 sym¬ 
metry. If this happens before the susceptibilities of the 
primary fields soften at Qj, then the Z2 rotational sym¬ 
metry breaks prior to other symmetry breakings. We 
emphasize that the nematic order parameter p involves 
the combination -|- 4)2, hence one cannot separate 
SDW induced and iCDW induced nematic order, even 
when Xs and Xc are not equivalent. 

We solve the action in ([3]) within the saddle-point ap¬ 
proximation, similarly to what was done in Refs. 0- We 
find that at 0s,Cc ~ 0, a non-zero nematic order pa¬ 
rameter {p) yf 0 emerges when the correlation length 
02 = ir/g, or, to logarithmic accuracy in g <C 1, at 
Ts = 27rps/| loggi, where ps is the stiffness of the 0(4) 
non-linear a model associated with Eq. ([2]) . It is instruc¬ 
tive to compare this result with the case where only 0(3) 
SDW fluctuations are present. In that case, the nematic 
order emerges when 30 q^3j = ^T^/gop,)^ and the transi¬ 
tion temperature is Tg = 27r/5s/| log,y 5 o( 3 )l, where gop,) 
is the coupling in the SDW 0(3) model. As a result, 
to obtain the same Tg, one needs a much smaller cou¬ 
pling constant 50 ( 3 ) ~ 5o(4)' Consequently, at T = Tg, 
the correlation length 0o(4) in the 0(4) case is propor¬ 
tional to 00 ( 4 ) a/ 0 o( 3 ), i-e. it is much smaller than it 

would be if nematicity was driven solely by SDW fluctu¬ 
ations. This is consistent with NMR [la and neutron 
scattering data [2l| in the paramagnetic phase of FeSe, 
which point to the presence of SDW fluctuations, albeit 
weaker than in the Fe-pnictide compounds. The rapid 
increase of the correlation lengths below Tg, obeying 
072 = jg g^jgQ consistent with the increase 

of 1 /TiT and the inelastic neutron signal (Tsl . 

The 0(4) Ising-nematic scenario also addresses why 
no magnetic order appears down to the lowest temper¬ 
atures. The SDW and iCDW orders compete yia the 
bi-quadratic term {u — g) M 2 4)2 in the low-energy action 
of Eq. ([2]). As a result, for U 2 > 0, fluctuations of the 
sub-leading iCDW channel suppress the transition tem¬ 
perature of the leading SDW channel. Such a suppression 
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FIG. 3: Density plot of the nematic transition Tnem as func¬ 
tion of the bare iCDW and SDW transitions Ticdw and Tsdw- 
To mimic the effect of pressure, they start at the same neg¬ 
ative value —6 at zero pressure, for which the nematic tran¬ 
sition temperature is Tnem.o, and then vary in opposite ways 
upon increasing pressure, Ticdw < —d and Tsdw > —d- 

is the largest when the difference between the coupling 
constants |rsdw —Ticdw| is the smallest, which happens 
when the system flows towards 0(4) symmetry within 
RG, i.e. when Ep is small, such as in FeSe. 

Experimental signatures. We know discuss the ex¬ 
perimental consequences of the Ising-nematic order. The 
breaking of the Z 2 symmetry between the j = 1 and 
j = 2 components of the 0(4) field implies the break¬ 
ing of O 4 lattice rotational symmetry down to O 2 . This 
instantaneously triggers structural order due to the cou¬ 
pling to lattice. To investigate how Z 2 order affect the 
electronic states, we return to the original four-pocket 
model (with fermions near the two hole pockets described 
by the operators ci^k and C 2 ,k) and include the explicit 
angle-dependence introduced by the matrix elements for 
the transformation between orbital and band basis. This 
transformation has the particularly simple form ci^k = 
dxz cos 0 k — dyz sin 0 k, C 2 ,k = dxz sin 0 k + dyz cos 0 k if one 
considers circular hole pockets and neglects the dxy or¬ 
bital component on the electron pockets [ 1 ^ 

The feedback effect of the Ising-nematic order on the 
fermions takes place via the self-energy corrections in¬ 
volving the unequal susceptibilities of the primary SDW 
and iCDW fields at momenta Qi and Q 2 . These 
corrections not only shift the chemical potentials of 
the /i and /2 electron pockets in opposite directions 
(/i,k/i,k) - (/ 2 ,k/ 2 ,k) oc (v?), but also give rise to a 
d-wave like distortion of the ci and C 2 hole pockets: 
(cLkCi.k) - ( 4 ,kC 2 ,k^ oc (v?)cos20k (see FigHb). In 
the orbital basis, the latter corresponds to ferro-orbital 
order {d^^dxz) — {d\jzdyz) oc {(p) [l^. Note that besides 
the changes in the dispersions proportional to {p), there 
is an overall shift of the chemical potential, symmetric 
for the two electron and the two hole pockets. 


The behavior of hole pockets in the Ising-nematic sce¬ 
nario is consistent with the existing ARPES data that 
show a d-wave type elongation of one of the hole pock¬ 
ets, whereas the other hole pocket sinks below the Fermi 
level [ 1 ^. The behavior of the electron pockets in the 
2-Fe Brillouin zone is also consistent with the splitting of 
the chemical potentials of the fi and /2 fermions. 

We also investigate how pressure affects the nematic 
transition temperature Tg. Within our approach, Tg is 
defined by the condition 3^^= 47r/g. Upon pressure, 
the Fermi pockets become bigger, and the Fermi energy 
increases. As a result iCDW becomes less competitive 
and decreases, while ^g increases. The combination of 
these two opposite tendencies in general gives rise to a 
non-monotonic behavior of Tg. This is illustrated in Fig. 
[3] using a simple modeling in which « T — Tj, with 
Tj denoting the bare transition temperatures for SDW 
and iCDW (see caption). 

Note that in our analysis so far we considered U 2 (0) > 
0. If on the other hand this interaction is attractive, 
^ 2 ( 0 ) < 0, the iCDW phase is the leading instability, 
and the ground state manifold is Z 2 x Z 2 . In this case, 
the nematic and iCDW transitions are expected to be 
simultaneous 0- Although at present no microscopic 
mechanism is known to give U 2 (0) < 0 Q, this could be 
another possibility to explain the existence of nematic 
order without magnetic order in FeSe. Such an iCDW 
phase could be detected via its time-reversal symmetry 
breaking, which would be manifested in, e.g., /iSR mea¬ 
surements. The phase diagram under pressure can be ex¬ 
plained by assuming that under pressure U 2 would change 
sign and SDW would become the leading instability. 

Summary In summary, we propose a natural exten¬ 
sion of the Ising-nematic scenario to explain the puzzling 
nematic state observed in FeSe. Our scenario relies on the 
smallness of Ep and explains the onset of nematic order 
far from magnetism due to the near degeneracy between 
the SDW channel and an iCDW charge-current density 
wave channel. This near-degeneracy could result in the 
nucleation of local iCDW order in the presence of point¬ 
like impurities, which favor iCDW against SDW order 
[ 4 ^ 1 . While these fluctuations cooperate with magnetic 
ones to break the tetragonal symmetry, they compete for 
long-range order and reduce both T/v and the magnetic 
correlation length at the onset of nematic order. We 
argue that this Ising-nematic scenario can also explain 
the observed non-monotonic dependence of the nematic 
transition temperature Tg upon pressure. 
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